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ABSTRACT 

The origin and sustenance of large scale galactic magnetic fields has been a long stand- 
ing and controversial astrophysical problem. Here an alternative to the "standard" a — O 
mean field dynamo and primordial theories is pursued. The steady supply of supernovae 
induced turbulence exponentiates the total field energy, providing a significant seed mean 
field that can be linearly stretched by shear. The observed micro- Gauss fields would be 
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produced primarily within one vertical diffusion time since it is only during this time that 

^ ■ 

^p" 1 . linear stretching can compete with diffusion. This approach does not invoke exponential 
O 

mean field dynamo growth from the helicity a-effect but does employ turbulent diffusion, 

a' 

which limits the number of large scale reversals. The approach could be of interest if the 



non-linear Galactic dynamo helicity effect is suppressed independently of the turbulent 
diffusion. This is an important, but presently unresolved issue. 

Subject Headings: Galaxy: general; galaxies: magnetic fields; ISM: magnetic fields 
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1. Introduction 



Magnetic fields are important to dynamics or emission in almost all astrophysical sys- 
tems. The formation of observed micro-Gauss large scale magnetic fields in the interstellar 
medium (ISM) of spiral galaxies has been considered a fundamental unsolved astrophysical 
problem (Beck et al. 1996; Zweibel & Heiles 1997). There is an important distinction be- 
tween the large and small scale magnetic fields observed in the ISM. Observations indicate 
that the small scale field is typically ordered on sub-kiloparsec scales and superimposed 
on that is a large scale, primarily azimuthal field. While it is more or less agreed that the 
small scale component of this field is injected and sustained by supernovae turbulence (e.g. 
Beck et al. 1996), the generation of the large scale field is where the controversy lies. 

Whether or not the large scale fields are of primordial origin or are produced in situ 
is difficult to determine observationally, and the "standard" competing explanations both 
suffer from difficulties. For example even if primordial fields could be produced, in situ 
processing of the field during the galactic lifetime needs to be addressed, and would likely 
dominate the present character of the observed field. The standard in situ mean field "a — O 
dynamo" (Parker 1979; Moffatt 1978; Ruzmaikin et al. 1988) has been threatened by some 
non-linear simulations (e.g. Tao et al. 1993; Cattaneo & Hughes 1996). While stretching 
of field lines by differential rotation (which provides the "^-effect") is not controversial, 
the helical property of the turbulence (which provides the "a-effect") and the turbulent 
diffusion (the " (3 effect") may be suppressed. Here I suppose that the a-effect is suppressed 
but that the (3 effect is not. Whether these two effects can be disentangled is unclear (Field 
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& Blackman 1997) but the possibility can still be explored. 

Turbulent amplification of the small scale field can provide a steadily supplied seed 
mean field (MF) whose subsequent stretching by differential rotation may produce and 
sustain large scale Galactic magnetic fields without a dynamo a-effect or a primordial seed. 
I first give the basic magnetic MF equations and summarize the complications of standard 
MF dynamo theory. The MF equation is then solved without the ct-eflect, but keeping in 
turbulent diffusion. It is shown that in 1 vertical diffusion time, differential shearing can 
increase the MF by an order of magnitude. It is also shown that the MF is maximized on 
a radial scale of a kpc or so: Turbulent diffusion tends to favor large scale field stretching, 
but this competes with the seed field's inverse dependence on scale. Simple statistical 
arguments are then used to predict the most likely number of reversals. That the small scale 
field is not necessarily dominant on only one scale is also emphasized. The assumptions 
employed here are no more controversial than in standard MF theory. Generally speaking, 
MF theory as applied to the Galaxy is most certainly an oversimplification, but it does 
provide a useful framework from which some understanding can be gained. 

2. Aspects of Standard Theory 

Writing the magnetic field (and velocity) as a sum of mean and fluctuating compo- 
nents, ie. B = B+B', the MF equation can be obtained by spatially averaging the magnetic 
induction equation over scales large compared to the fluctuations but small compared to 
the overall scale of the system (e.g. galaxy): 

9 ( B = Vx(vxB)+Vx(v'xB')+i/ M V 2 B, (1) 
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where v is the velocity and vm is the magnetic diffusivity, and the primed and barred 
(or bracketed) quantities indicate fluctuating and mean values respectively. The turbulent 
EMF in (1) can be written (v' x B') = a ij (B,v)B j - /3 ijk (B,v)VjB k + - where 
indicates higher order gradients. In linear kinematic dynamo theory for isotropic incom- 
pressible plasmas, a oc (v' • / V x v'(t')dt'} and oc (V • / v'{t')dt'). In the dynamic, 
non-linear theory, these can be functions of B, and are not necessarily isotropic. 

In the standard a — O dynamo (aQD) applied to the Galaxy (e.g. Ruzmaikin et al. 
1988), supernovae (SN) induced turbulent eddies on ~ lOOpc scales stretch field lines into 
loops or cells. The coriolis force, in principle, conspires to statistically twist all of these 
loops in the same direction, providing a much larger scale mean loop. This, resulting 
from the large scale reflection asymmetry, provides the non-vanishing pseudoscalar a- 
effect (Parker 1979; Moffatt 1978). The outer portions of these loops must incur turbulent 
diffusion to leave a net mean flux in each hemisphere of the Galactic disk (Parker 1979). 
The large scale field formed in this way is further sheared by differential rotation (the 
fi-effect), providing new toroidal field and starting the process again. In principle, this 
feedback leads to exponential MF growth of a primarily azimuthal field with a growth time 
~ 4 x 10 8 yr (Ruzmaikin et al. 1988). 

Standard kinematic aVtD treatments ignore the backreaction of the growing magnetic 
field on the turbulence. Because high magnetic Reynolds numbers make the last term in 
(1) negligible on the energy containing scales, the field exponentially grows to equipartition 
with the turbulent energy by a fast dynamo (FD) on a time scale much shorter than any 
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MF evolution time (e.g. Parker 1979) and does not require helicity. In combination with 
even a weak MF that is > p 1 / 2 v'/Rl( 2 (where R m is the magnetic Reynolds number, and p 
is the density), this may make Lorentz forces (Kulsrud & Anderson 1992; Cattaneo 1994) 
lock a significant fraction of motions into oscillations. The magnetic fields act like springs 
that the required turbulent motions must fight against. Although simulations in 2-D show 
P suppression (Cattaneo 1994), there are not yet simulations that show (3 suppression in 
3-D. There have been some simulations showing a-effect suppression in 3-D (Tao et al., 
1993; Cattaneo & Hughes 1996), and others that do not (Brandenburg & Donner 1997). 
Intermittency (Blackman 1996; Subramanian 1997) and the nature of the forcing function 
may play a role in overcoming both a and (3 suppression in 3-D. Basically, there is no clear 
consensus on what happens in fully non-linear mean field dynamo theory with respect to 
the backreaction, even as to whether suppression of a and (3 are intertwined (e.g. Field & 
Blackman 1997) or independent. 

3. Linear Mean Field Growth from a Random Seed Field 

SN inject turbulent energy into the ISM, and also inject magnetic field (Ruzmaikin et al. 
1988; Rees 1994). Because of the observed dispersal of heavy elements (Rana 1991), the 
supernova ejecta at least mix with the remnant material. Theoretical estimates for the 
mean seed field injected from SN range from 10~ 13 Gauss from simple flux freezing, to 
10~ 8 Gauss for including winding in pulsar winds (Rees 1994). 

Both theory and simulation (Parker 1979; Piddington 1992; Beck et al. 1996) show 
that the FD builds up small scale magnetic field energy on a growth time of order the energy 
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containing eddy turnover time ~ l/v' ~ 10 7 yr, where / ~ lOOpc is usually taken as the 
energy containing eddy scale and v' > lOkm/sec is a typical observed speed of these eddies. 
Thus, once the Galactic volume is full of a weak seed magnetic field from the first set of 
SN, the next generation of eddies stirs the field to equipartition. The SN remnants fill the 
Galactic disk (height ~ 500pc by radius 12kpc) every 10 7 years given their observed rate of 
~ 0.02yr _1 (Ruzmaikin et al. 1988), and this maintains a steady random field energy. How 
the field actually mixes from the SN to the ambient ISM is complicated. The amount of 
magnetic annihilation, the amount of enhancement, the geometry /topology of the injected 
field (e.g. Ruzmaikin et al. 1988), and the role of boundary instabilities are all subtle issues. 
Despite these complications, the basic picture of SN seed field injection and subsequent 
stirring as described above, leading to equipartition fields at a time < 2.5 x 10 8 yr (where 
this upper limit comes from using the 10 _13 G seed value given above) and all subsequent 
times, is consistent with observations: The magnitude of the random field is observed to 
be B ~ 5 - 10 x 10" 6 G (e.g. Heiles 1994; Rand & Kulkarni 1989). Rand & Kulkarni 
(1989) impose a single cell model whose best fit small scale over which the field is ordered 
is 50-100pc. Though this has become the standard quoted range for the small scale, it will 
be emphasized later why multiple and larger cell sizes (Ferriere 1994) are important. 

The above small scale field would give a corresponding MF of magnitude B ~ B/N 1 / 2 
where N is the number of small scale coherence volumes in the region of averaging. Such 
a residual large scale field has been argued to be a viable source of seed field for the 
Galactic aQD (Ruzmaikin et al. 1988; Rees 1994), but below I suggest that even if a is 
suppressed, an appropriate large scale field can still be produced. Assume a is suppressed 
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well below the critical value required (Parker 1979) for the standard a — O dynamo growth, 
ie., a « a cr it ~ /3 2 /(fi/i 3 ) (where h is the disk height), so we can then ignore it in what 
follows. The MF induction equations in cylindrical coordinates become 



d t B r = (3V 2 B r 



(2) 



dtB^ = rd r QB r + f3V 2 B (j}: 



(3) 



d t B z = /3V 2 B 



(4) 



Eqs. (2) and (4) are decoupled from (3), implying pure diffusion of B r and B z . As 
in standard treatments (e.g. Ruzmaikin et al. 1988), the O-effect (differential rotation) 
increases the azimuthal field in (3) linearly on a time of order the rotation time. For the 
Galaxy this is O -1 ~ 3.3 x 10 7 yr. Because the rotational energy far exceeds that which 
can be transferred into the magnetic field during the age of the universe, the 0-effect is 
not controversial; there is no back-reaction on the large scale rotational motion. As shown 
below, this O-effect can generate a factor of ~ 10 increase in the large scale field. This is 
sufficient without an a-effect since the seed field is continually supplied. 

To solve the equations and determine the dominant scale of the MF, I assume that 
B = B t e ( -* k ' x - > , where k is the wave vector of the MF and the subscript t labels the time 
dependence. The MF equations (e.g. Ruzmaikin et al. 1988) without an a-effect and with 
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homogeneous (3 for the azimuthal and radial fields are then 



d t B (j)t = B rt fQ-(3k 2 B cj)U 



(5) 



and 



d t B rt = -/3k 2 B rU 



(6) 



where fQ = rdVL/dr. Solving (6) gives B rt — S r oExp[— k 2 



fit], so (5) gives 



dtB+t = fnB r0 E w [-k 2 /3t] - k 2 (3B rt . 



(7) 



Multiplying both sides of (7) by Exp[/c 2 /3t], using the chain rule and solving gives 



where I have taken B r0 ~ B^ ~ Bo/3 1 / 2 = B/(3iV) 1/2 , N ~ hDR/(4nl 3 /3), and h is 
the Galactic scale height, while D and R are the azimuthal and radial mean field gradient 
lengths corresponding to their wave vectors, and defined only for scales > h. 

The dominant contribution to the MF at any one time is that produced within a 
vertical diffusion time, r^, from the observation time. Thus fit ~ l//i 2 , and from (8) 



B<t>Td V = B(fnh 2 /[3 + l)[3/i J Di?/(47r/ 3 /3)]- 1 /2 Ex p[-(l + h 2 / D 2 + h 2 /R 2 )]. (9) 



The scale height of the disk is fixed at h ~ 500pc, but B^ Tdv can be extremized as a 
function of R and D, giving a maximum at R = D = 2(f3t)~ 1 / 2 . Using f3 ~ 10 26 cm 2 /sec 



B^ t = B(fQt + l)(3iV)- 1 / 2 Exp[-(A; 2 + k 2 z + k 2 )(3t] 



B(fQt + l)[(3hDR)(4TTl 3 /3)]- 1/2 Exp[-(D- 2 + h~ 2 + R~ 2 )/3t], 



(8) 
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(e.g. Ruzmaikin et al. 1988) a Galactic disk scale-height of h ~ 500pc, ~ 10 15 sec 1 , 
and / ~ 1 this gives D = R = lkpc. Thus B^ Tdv = 1.42 x 10" 6 (B/5 x 10" 6 G)G . If instead 
we take D = 3kpc, to match the radial scale measured by Faraday rotation (e.g. Rand & 
Lyne 1994), this becomes 1.02 x 10~ 6 (S/5 x 10~ 6 G)G. The approximate magnitude of the 
local MF (e.g. Heiles 1994) may therefore be reproduced without the cn-effect. 

Note that B$ Tdv depends on the value of the characteristic averaging azimuthal dis- 
tance as D _1 / 2 Exp[— h 2 /D 2 ], on the characteristic small scale structure size / to the —3/2 
power, and linearly on B. The actual small scale field of the Galaxy has been shown from 
observations to be inconsistent with a single scale size (Rand & Kulkarni 1989): The sta- 
tistical dispersion between the observed field and their model large scale field shows little 
evidence of fall off with the distance to pulsars as it should if the the single cell size model 
were appropriate. This highlights the importance of super-bubbles and other larger scale 
fluctuations known to be important to the field structure (Ferriere, 1996). If the region 
over which the measured MF were composed of primarily > 200pc instead of lOOpc struc- 
tures, then the estimate given above is magnified by an additional factor > 2 3 / 2 = 2.8. 
Also, Heiles (1994) finds that the total magnetic energy does not scale simply with the 
mean azimuthal field as measured by Faraday rotation in different parts of the Galaxy. 
This can be explained in the present model, since the mean field is proportional to the 
RMS field divided by iV 1 / 2 , where N is the number of small scale cells of uniform field in 
the region determining the Faraday rotation measure. Regions of different cell sizes would 
therefore produce different observed mean fields even if the total magnetic energy density 
were the same. All of this highlights the possible importance of multiple small scale sizes. 
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4. Discussion of Mean Field Reversals 



The calculation of section 3 shows that a scale height of h = 500pc, maximizes the az- 
imuthal MF for a radial averaging scale of D = R = lkpc. This results from two competing 
effects: [1] The time for shear to increase the field strength by an order of magnitude is 
relatively independent of scale. Thus large scales are preferentially sheared in a fixed time 
because the competing turbulent diffusion depends on scale squared. [2] However, the 
initial seed field depends inversely on the averaging scale. The scale of lkpc optimizes 
[1] and [2]. This defines the minimum radial scale over which the maximum average az- 
imuthal field could reverse sign. This does not mean that there would be necessarily be 
reversals every lkpc. It means that between lkpc annuli, the mean field may or may not 
reverse. Within a lOkpc Galactic radius, there are ~ 9 interfaces between lkpc annuli. 
The probability P(n) of observing n reversals by Faraday rotation would be 9 'choose' n, 
i.e. 9!/(9-n)!n!, which is maximized for n=4 or 5. 

Galactic Faraday rotation observations can determine the sign of the large scale field 
in the line of sight (cf. Beck et al. 1996; Zwiebel & Heiles 1997). (Unlike Galactic 
measurements, where pulsar dispersion measures can be used, extragalactic measurements 
require independent determinations of the density to obtain any information from Faraday 
rotation. The data for external galaxies are therefore less reliable (Heiles 1994; Zwiebel 
& Heiles 1997).) Generally, a large scale theoretical Galactic field model is statistically 
compared to observations (e.g. Rand & Kulkarni 1989). Field reversals seem to occur in 
each of the two interarm regions immediately inside of the solar circle (Beck et al. 1996; 
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Heiles 1994) with perhaps two more outside. The reversals are not necessarily periodic 
between spiral arms (Vallee 1996). Also, because of fluctuations in rotation measure data 
for some quadrants (Rand & Lyne 1994; Beck et al. 1996) averaging over smaller scales 
then shows smaller intermediate scale reversals. This again highlights that intermediate 
scales (Rand & Kulkarni 1989) from 50-500 pc complicate theoretical and observational 
interpretations. The precise structure of the large scale field in spiral galaxies is difficult 
to conclusively determine (Beck et al. 1996). 

Note that turbulent diffusion is distinct from dissipation. The former describes a 
transfer of magnetic energy between scales, whereas dissipation is a removal of magnetic 
energy. Turbulent motions on sub-kiloparsec scales both randomize the mean field and 
amplify the small scale field, thereby re-seeding the mean field. Though a turbulent cascade 
drains energy to the dissipation scale, the magnetic energy is steadily replenished by the 
FD and the total magnetic energy density remains steady. 

Previous work has recognized the importance of diffusion for reversal reduction (Poezd 
et al. 1993). In fact, the weaker the a-effect in dynamo models, the less vigorously the 
aVtD can compete with turbulent diffusion and the fewer reversals that survive. Primor- 
dial models are sometimes employed with the assumption that turbulent diffusion is not 
operating (Zweibel & Heiles 1997). Though this seems unlikely, other proposed mecha- 
nisms would then be needed to eliminate reversals (c.f. Zweibel & Heiles 1997). Another 
possibility is that the winding of a proto-galactic field in the subsequently formed galaxy 
(Howard & Kulsrud 1997) generates the correct number of reversals. 
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For some external galaxies (e.g. NGC6946) , observations indicate that the large scale 
field is actually stronger in the interarm regions (Beck & Hoernes 1996). In the present 
approach, the deficit of large scale field in the spiral arms would be the result of a reduced 
shear there (Elmgreen, 1994) and thus an / < 1 in Eq. (5). This is generally consistent with 
rotation curves of NGC6946 and other galaxies which show reduced differential rotation in 
spiral arms (Sofue 1986; Rubin et al. 1980). In contrast, an enhanced MF strength might 
result in the arms if their increased electron density dominates the effect of reduced shear. 
The total magnetic energy can be larger in spiral arms if the turbulent energy is higher 
there. Varying density complicates the interpretation of rotation measures of external 
galaxies if the density variation cannot be independently measured. 

5. Discussion 

It is important to understand whether a and (3 can actually be disentangled. If so, the 
main point herein is to suggest that it may not be absolutely certain that the observed 
large scale Galactic magnetic field requires a dynamo a-effect even if the mean field is 
produced in situ. The well-known exponential growth of small scale field by the FD and its 
steady replenishing of seed MF, combined with the subsequent linear growth of large scale 
azimuthal field by the fi-effect might supply a large scale Galactic field without requiring 
the dynamo a-effect. The linear growth may be sufficient because it proceeds faster than 
the time for the field to diffuse below micro-Gauss values. The observed field of any spiral 
galaxy would be that produced within Td v of the time of observation, and since the field is 
steadily replenished, this statement is true at any time in a galaxy's lifetime > 10 8 yr from 
the time of the galaxy's birth. The most likely number of reversals in the large scale field 
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within lOkpc radius would be of order 4-5 in the presence of turbulent diffusion for the 
simplest approach. (Unlike the cellular model of Michel & Yahil (1973), here turbulence 
is important.) If flux tubes were present with a small volume filling fraction, and/or if the 
the energy containing small scale of the field were much reduced from the semi-empirically 
determined lOOpc scale, too many mean-field reversals might be produced by the present 
approach. A small filling fraction may also aid the standard dynamo (Blackman 1996; 
Subramanian 1997), making the approach herein less useful. However, unlike the Sun, 
the tube filling fraction in galaxies may be large if the average particle pressure does not 
overwhelm the magnetic pressure (Blackman 1996). 
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